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S ' Abstract 

The 7T° — ► 77 decay width is analyzed within the combined framework of Chiral Perturbation 
Theory and the 1/N C expansion up to 0(p G ) and C(p 4 x 1/N C ) in the decay amplitude. The rf 
is explicitly included in the analysis. It is found that the decay width is enhanced by about 4.5% 



due to the isospin-breaking induced mixing of the pure U (3) states. This effect, which is of leading 



order in the low energy expansion, is shown to persist nearly unchanged at next to leading order. 
The chief prediction with its estimated uncertainty is r„-o^ 77 = 8.10 db 0.08 eV. This prediction at 
the 1% level makes the upcomming precision measurement of the decay width even more urgent. 
Observations on the 77 and rf can also be made, especially about their mixing, which is shown to 
be significantly affected by next to leading order corrections. 

PACS numbers: 11.30.Rd, ll.40.-q, 12.39.Pe, 13.20.Cz 



I. INTRODUCTION 



In the chiral SU(2) limit (m u> d = 0) the 7r° — > 77 decay amplitude is precisely known 
to order a[ttj. The amplitude is 0(p 4 ) in the low energy chiral counting. It is determined 
entirely by the anomaly induced on the divergence of the axial current A 3 ^ = q^f/y^q by the 
electromagnetic interaction, and is expressed in terms of the only two available quantities — 
the fine structure constant a and the pion decay constant F n — the decay width in this limit 
is thus given by r„-o^ 77 = (-p-) 2 (-^-) 3 - The explicit breaking of chiral SUl(2) x SUr(2) 
symmetry induced by non-zero u- and d-quark masses generates corrections to the chiral 
limit result, and it is the purpose of the present work to evaluate these corrections as well 
as to understand their origin. In order to achieve this goal it it is crucial to perform the 
analysis in the extended framework of three flavors supplemented by the 1/N C expansion 
in order to include explicitly the rf degree of freedom. In this three flavor framework the 
corrections turn out to be of two types: 

i) those due to isospin breaking (mixing corrections) that are proportional to (m u — 
m d)/ m s or to N c {m u — m d )/A x , both giving contributions to the decay amplitude 
that, according to the counting defined in the next section, are 0(p 4 ), i.e., of the same 
order as the leading term, and 

ii) those proportional to m u ^/A x that are of subleading order — 0(p 6 ) — and which stem 
from different sources, as shown below. 

The inclusion of such corrections is crucial for a prediction of the 7r° — ► 77 width at the 
1% level, which is the level of theoretical precision required by the forthcoming dramatic 
improvement expected in the experimental measurement of the tt° width via the Primakoff 
effect. The PRIMEX experiment at Jefferson Lab |2j is aiming at a measurement with an 
error about 1.5%, which is several times smaller than the 7.1% uncertainty in the current 
world-average value 7.74 ± 0.55 eV ||. However, this quoted experimental uncertainty is 
open to question as can be seen by the large dispersion of the experimental results which 
suggests that the errors of the individual experiments have been underestimated. Indeed, 
a direct measurement gives T^o = 7.25 ± 0.23 eV||], a production experiment in e + — 
collisions yields T^o = 7.74 ± 0.66 eV||, while Primakoff effect experiments, all dating back 
to the earlv seventies, give disparate values: a large number — 11.61 ± 0.55 eV HI, and two 



that are consistent with the current world average — 7. 22 ±0.55 eV [|7| and 7.93 ±0.39 eV ||. 
This unsatisfactory experimental situation together with the rather precise theoretical pre- 
diction derived in the present work clearly lend great significance to the upcoming PRIMEX 
measurement. 

Within the two-flavor framework, wherein the strange quark is integrated out, any correc- 
tions to the amplitude are 0(p 6 ) and reside entirely in F n o |J or in the 0(p 6 ) odd-intrinsic 
parity Chiral Lagrangian||l 1||, also known as the 0(p G ) Wess-Zumino (WZ) Lagrangian. At 
leading order, the ensuing theoretical prediction (taking F n o = F n +) is r 7r o_ >77 = 7.725 eV, 
a result that agrees well with the experimental world average within its generous error. 
The analysis within SU(2), however, does not provide insight on the origin of the 0(p 6 ) 
WZ contribution just mentioned. Such an insight can be gained by instead performing the 
analysis in the t/iree-flavor framework, as has been shown by Moussallam[|l^]. In particu- 
lar, he pointed out that the primary corrections to the tt° width result from the leading 
order isospin breaking effects mentioned above in i), which stem from the m u ^ m^-induced 
mixing between the pure isospin state ir° and pure SU(3) states rj and rf (to be denoted 
below by 7C3, irg and ttq respectively). The present work confirms that such isospin breaking 
corrections persist as the dominant effect when next to leading order (NLO) corrections are 
included. 

In this work then the tt° — > 77 decay rate is evaluated to NLO within £/l(3) x Ur(3) 
chiral perturbation theory wherein the 7/ meson is included consistently by means of the 
1/N C expansion since in the large N c limit the rf becomes a Goldstone Boson. Such a 



framework was recently developed by Herrera-Siklody, Latorre, Pascual and TaronfT3|| and 
by Kaiser and Leutwyler fT4|, |T5f , who showed that a simultaneous chiral and 1/N C expansion 
leads to an effective theory for the pseudoscalar nonet that is not only internally consistent 
but is also very useful in practice, as the present work shows. 

The chief result of our paper is that the 7r° — > 77 width is enhanced by about 4.5% from 
the lowest order chiral anomaly prediction, a result expected to hold within an uncertainty 
of ± 1% after NLO contributions are included. The magnitude of this enhancement agrees 
with that obtained in the analysis of Moussallam [12], where the NLO corrections were not 
implemented in a consistent fashion, as done in the present work. 



II. TWO-PHOTON DECAY AMPLITUDES 



The decay amplitudes of 7r°, rj and rf into two photons can be obtained from the Ward 
identities satisfied by the three axial vector currents A ^ = 5(77^75 X a q (a = 3, 8, 0), where A a 
are U(3) generators (A being the U(l) generator) normalized via (X a X b ) = Tr(X a X b ) = 25 ab . 
In the presence of the strong and electromagnetic interactions, the divergence of the axial 
vector current is given by: 

Q * A l = -J^^ a Q 2 }FF + ^-(X a }GG+-q l5 {X a ,M q }q + --- } (1) 

where A4 q is the quark mass matrix and eQ is the electric charge operator. Here 
FF = -jepvpaF^F 9 ® , F^ u being the electromagnetic field tensor, and similarly GG = 
\t l iu P aG atJ ' u G apa , G^y being the gluon field [1] , and the ellipsis denotes terms irrelevant 
to this work. 

The two-photon amplitudes can be obtained by considering the matrix elements 

( 77 I d»Al I 0) = C a ^<77 I FF | 0) + 5 a0 ^<77 I GG | 0) 

+ l(ll\qid^ a ,M g }q\0), (2) 



where C 3 = 1, C 8 = l/y/3, and C = y/8/3. 

If p denotes the total momentum of the final two-photon state, then in the limit of small 
p 2 Eq. H] admits a low energy expansion and can be expressed as: 

X>7 I ^,p)(^P I pai I o) p2 _^ |+ . £ = Ca^(m I FF 1 0) 
+ *«o^E<ry I «*p)fa, P I gg 1 o> p2 _ j ^ + . e 

+ \ X^ 77 I n a,P)(^a,P I qih{X a ,M q }q I 0) 2 fl l 2 ,,- +••• ' 

z p iW„ -+- It 

where (77 | Tr a ,p) are the two-photon amplitudes, the ellipsis denotes contributions from 
excited mesons as well as from the continuum, all being of NLO or higher, and the mass 
eigenstates n a that correspond to the physical 7r°, 77 and rf are given by 

7T« = ^ Qaa^a , (4) 
a=3,8,0 



Ml r Tlirmio'Vimit tnr> pnnvpTitinriQ in RinrVnn flnH T^mll ftrn iiqpH 



where the mixing matrix that diagonalizes the mass matrix is parametrized in terms of Euler 
angles 6*3, 6*8 and —9 



^ C 3 C 8 — Cq s 3 s S C3S8 + C8C0S3 — S3S0 ' 







— C 8 S 3 — C3C0S8 ~ S 3 S 8 + C3C8C0 — C3S0 

y -s 8 s c 8 s c j 

where q = cos^j and Sj = sin#j. Here, for small mixing, the projection of the physical 
7T° onto 7r 8 is given by the angle e ~ 9 3 + 9 8 , the projection of the physical 77 onto n 
is approximately given by — 6q (8 can therefore be identified with the well known 77 — rf 
mixing angle as it is customarily defined), and the projection of the physical 7r° onto 7Tq is 
given by e ~ —O^Oq. 

The NLO — i.e., 0(p 6 ) — corrections in Eq. [| reside in the terms displayed explicitly 
through their dependence on the masses and decay constants, as well as in pieces that stem 



from continuum and excited states. In works that preceded that of Moussallam [12], such 



as references II and flQfl , such mixing corrections as well as the NLO effects of the latter 
kind were disregarded. Ignoring such effects implies that only NLO corrections which are 
absorbed into the 7r° decay constant remain, as it was shown in |J; in that case and taking 
F^o = F n +, the predicted 7r° width is the previously mentioned 7.725 eV. As shown here, 
however, disregarding mixing in particular constitutes a very poor approximation. In the 
presence of mixing the pseudoscalar decay constants form a 3 x 3 matrix defined by the 
matrix elements of the axial vector currents, which connect the pseudoscalar mesons to the 
vacuum: 

<7r a ,p|^|0) = -z^F ao . (5) 

Indeed this decay-constant matrix contains all that is needed to calculate the three two- 
photon amplitudes, except the contributions stemming from the 0(p A ) WZ Lagrangian, and 
its evaluation is the centerpiece of the present work. 

As already mentioned in the introduction, the present analysis includes the 7/ as an ex- 
plicit degree of freedom, which, in order to be consistently implemented in an effective theory, 
requires the validity of the 1/N C expansion, wherein, taking the chiral SUl(3) x SUr(3) limit, 
= 0(1/N C ). Thus, in the framework of the 1/N e expansion M3 should be considered as 
"small," and its explicit inclusion becomes consistent with having a simultaneous low energy 



chiral expansion. This explicit inclusion of the rf in the low energy expansion implies that 
M3 must count as a quantity of 0(p 2 ), which in turn implies that l/N c should be counted 
as a quantity of the same order. Indeed, a consistent effective theory can be formulated with 
such a counting scheme [O, IT3, IH3], and it is interesting to note that 1/N C and the magnitude 



of SU(3) breaking are comparable in size in the real world. 

Taking the chiral limit — M. q — > — in Eq. ||] and neglecting the electromagnetic piece — 
i.e., the FF term — , equating the residues at p 2 = M 2 , lead to the well known relation 

Ml = ^£jrW\GG | 0), (6) 

with M being the rj' mass in the chiral limit. Here the lowest order result 

(ir- a ,p\d»A« \0) = 5 aaP 2 F 

was used, where Fq is the pion decay constant in the chiral limit. In the large N c limit, Fq 
scales as \J~N~ C , while in the chiral limit the ratio F /F v ' is equal to unity up to corrections 
of order 1/N C . 

On the other hand, for non- vanishing quark masses equating the residues in Eq. [5] yields 

<*a,P | qi^\M q }q | 0) = -2i{{it- a ,p \ PA} | 0) 

- 5 a0 VQ^(7r- a ,p\GG\0)) (7) 

where p 2 = M 2 . As is well known, the p 2 dependence of the LHS appears first at 0(p 6 ) flIS| , 
which would affect the two-photon amplitudes at 0(p 8 ), i.e., beyond the accuracy needed 
in this work. Thus, it is consistent to use Eq. [7] at p 2 = M| to represent the LHS in the 
entire low p 2 domain. 

At LO — 0(p A ) — then, Eqns. |^, and|^ yield immediately the result for the two-photon 
amplitudes 

(77 I 7T S ) = "^^^(77 I FF | 0) (8) 

a=3,8,0 

At this order the decay constant matrix is simply given by 

Faa = QaaFo (9) 

where Sa is the mixing matrix obtained from the 0(p 2 ) mass formulas. Of course, the 

rpsnlt. nf P,n R pninpiHpR with thp rpcnilt nhtninpH hv mpuns nf thp (1(rfi\ W7 tprm inpliiHinc 



explicitly the singlet pseudoscalar. The purpose of carrying out the above Ward identity 
analysis is, however, to make more transparent the origin and structure of the higher order 
corrections and to set the stage for the inclusion of NLO modifications. 



III. LEADING ORDER RESULTS 

The leading order mass formulas are obtained from the 0(p 2 ) chiral Lagrangian. The 
U (3) field is parametrized by the unitary matrix: 

where F = 92.42 MeV at LO. The 0(p 2 ) Lagrangian with the standard definitions of 
covariant derivatives and sources x EH ^ s gi ven by: 



£ (2) = \f 2 (d,ud^} + \f 2 ( x u^ + x ] u) - \m^I (11) 

and the mass matrix in the tt 3 , tt 8 , tt sector of interest resulting from is: 



M 2 LO = B 



( 2m -^(m u -m d ) ^J\(m u -m d ^ \ 

^(m u -m d ) §(m + 2m s ) 

f2 



(12) 



\Jl(m u -m d ) -^{m s -m) ^ + §(2m + m s ). J 

Using the leading order mass formulas — e.g., M 2 + = 2B m, with 2m = m u + m d — and 
extracting isospin breaking from the K + — K° mass difference via Dashen's theorem to 
eliminate the EM contributions [3] , a best fit to the masses yields a singlet mass Mq of 
approximately 850 MeV, and the Euler angles: 8% = 1.57°, #8 = —0.56°, and 8q = —18.6°. 
This fit yields then e = 9% + 6$ ~ 1°, which is substantially larger than the value [|l6| 

a/3 m d - m u 



0.56° (13) 



that arises in the limit M — > oo, when only octet degrees of freedom are included. The 
LO mass matrix, however, gives a poor result for the masses. In particular the r] mass is 
too low by almost 50 MeV, a problem that is generic at LO in the low energy and 1/N C 
expansions Jl7 . 
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Using the relations 

M 3 

(77 I 7Ta> = «,(77 I FF I 0), T- a =| K- a | 2 (14) 

which connect the decay amplitudes and associated widths, the fitted parameters at LO lead 
to a 7t° — > 77 decay width of 8.08 eV, which is 4.5% larger than the leading order result 
r^o^^ = 7.725 eV obtained in the two-flavor framework wherein mixing effects are moved 
to NLO. It should be noted, however, that the two-photon widths of the 77 and rj' predicted 
in this leading order fit are too large, the first being 22% and the second 20% larger than 
the corresponding experimental values. One of the chief reasons for this disagreement is 
that SU(3) breaking in the pseudoscalar decay constants is not included at LO. The rj — rj' 
mixing angle 9 turns out to be ~ —18.6° at LO, and will be reduced by almost a factor of 
two when NLO corrections are included. The magnitude of the observed LO enhancement of 
the 7r° width is in line with the ratio of isospin breaking versus m s {{md — m u )/m s ~ 2.3%) 
and versus Mo (Bo(rrid — i7i u )/Mq ~ 1.1%), and is therefore not surprising. It is important 
to note that the corrections due to mixing with 77 and with rj are of the same sign and of 
similar magnitude. This point is the primary reason why the 7/ must be explicitly included 
for a full understanding of the mixing effects. Although 7/ — 7r° mixing is smaller than the 
77 — 7T° mixing, the bare singlet state has an intrinsic two-photon amplitude larger by a 
factor a/8 that compensates for the smaller mixing. Overall, however, the LO fit is poor, 
and dramatic improvement results when the NLO corrections are included, as shown in the 
following section. 



IV. NEXT TO LEADING ORDER ANALYSIS AND RESULTS 

At NLO the amplitudes receive corrections of two types — those that affect the decay 
constants and mixing angles, and those that stem from the presence of excited states and 
which are included in the 0(p 6 ) WZ Lagrangian. 

The first type of correction requires the determination of the masses and decay constants 
to NLO and can be obtained in the standard fashion by calculating the two-point functions 
of axial vector currents, where the relevant diagrams are shown in Fig. 1. Up to 0{p 2 ) and 
O(p /N c ) such two-point functions require only the effective Lagrangians and C^. Chi- 
ral loop corrections are 0(p 2 /N c ) = 0(p 4 ), and therefore beyond the precision of the present 



calculation, but such loop corrections will be calculated merely as a means to estimate the 
size of possible contributions from terms of that order and also as a test on the practical 
validity of the 1/N C expansion. The lowest order Lagrangian has already been given in Eq. 
|TT| , while the next to leading order Lagrangian has the form given by [|T^, *H, [T5|, [H| 

£( 4 ) = ... + L 4 (D^D»U)(x j U + rfx)+L 5 (D fl rfD»U(x j U + rf X )) 
+ L 6 ( X *U + U^ X ) 2 + LrtfU - U^ X ) 2 + Lsixrfxtf + h.c.) 
+ ^D^D^-^-MxU^-x^U) 

+ tL 18 V6D^ (D^x ~ D»U X *) + iL 25 V6ir (U^xrfx ~ U X ] U X ] ) + • • • , (15) 

where only the terms relevant to this work are included. In the presence of the SU(3) singlet 
axial vector source field a°, D^tiq = <9 M 7ro — Fq a° . At this point it is important to note that 



the singlet axial current has non- vanishing anomalous dimension |18fl , which implies that 
some low energy constants (LECs) as well as the singlet field itq must be renormalized (the 
corresponding renormalized quantities will thus depend on the QCD renormalization scale 



I^qcd) [14, 15 1 . Since the renormalization of the axial current is subleading in 1/N C , such 



dependence appears first at the level of £ (4) LECs. It has been found in particular that the 



LECs Ai, A2, Lis and L25, all of which are subleading in 1/N C , must be renormalized [14 
implying that the values of these LECs will depend on the value of the scale Hqcd- Other 
quantities such as F v > and the singlet 7r field must be renormalized as well and depend on 
Hqcd through the renormalization factor Za associated with the singlet axial current. It 
is very convenient to make use of the asymptotic freedom of QCD to set Hqcd arbitrarily 
large and give the values of the LECs in that limit. Indeed, the renormalization factor of 
the axial current Za evolves to a fixed point that can be taken to be Za = 1 as \iqcd 00 • 
All quantities given in the following that have a dependence on \xqcd are then taken in 
this limit. It is well known that the low energy constants L 5 and L s are 0(N C ), while L 4 
and L 6 are subleading and 0(N®) Jl6| and are needed in order to renormalize the one-loop 
contributions from C^ 2 \ With 7/ as an explicit degree of freedom, L 7 is also subleading 
in 1/N C . The renormalized pieces of the subleading LECs are therefore set to zero at the 
chosen chiral renormalization scale \i in our analysis. On the other hand, the LECs Ai and 
A2 are 0(1/ N c ) and the corresponding terms in the Lagrangian are 0(p 2 /N c ) = 0(p A ) and 
thus must be included in the calculation. The low energy constant Ai provides an 0(1/N C ) 

rnrrprtinn tn ttip n Hppfiv nnnst.ant wtiilp Kntti A, anrl affppt pntrips in thp mass matriY 



FIG. 1: Two-point function of axial currents to NLO. The last three diagrams involve the coun- 
tertem insertions from 

£(4) 



involving the rf at order 0(p 2 /N c ). Finally, the terms involving the LECs Lis and L25 
are of 0(p A /N c ) which is beyond the order of the present calculation and therefore their 
renormalized pieces are set to vanish as well. The renormalized LECs are defined in the 



usual MS renormalization scheme 116 



A ("» = £(^-^ 4 ' +1+r 'w») < 16 > 

The /3-functions Tj associated with the LECs Li, and A» associated with A, that result 



from the chiral one- loop calculation are given by |L3], |14|, [y|: r 4 = 1/8, T 5 = 3/8, T 6 = 1/16, 
r 7 = 0, T 8 = 3/16, T 18 = -1/4, r 25 = 0,A 1 = -l/S, and A 2 = 3/8. 

The two-point functions of axial currents can be written in momentum space as 

f d 4 x e*"<0 I T(A°(x)A b u (0) | 0) = p,p v £ Fg{p 2 )A- a {p 2 )F- aa {p 2 ) + •• ■ , (IT) 

where the term explicitly shown contains the light pseudoscalar poles and A a (p 2 ) is the 
propagator of the corresponding mass eigenstate. From the location of the low energy poles 
and the residues of the two-point function the light pseudoscalar decay constants and masses 
are extracted, yielding 



M 2 ab = M 2 LO ab - (a CT + oioop - Im 2 lo , ct + ^0 \) ab , a, b = 3, 8, 



M 2 a =M 2 LOa - {{a CT + o loov ) a - —M 2 LO a (<f) CT + -0) a ), a ^ 3, 8, (18) 

-^0 * 



The corresponding decay constants are given by: 



Faa ®abFba 

1 3 

Fab = F ba = F (5 ab - -((per + ^<P)ab), a, b = 3, 8, 
F a = Fo-h^cr + ^a, a + 3, 8, (19) 



where the following definitions were made: 



hcTab = -(^({A a ,A 6 }^} + AiW M ) 

Kb = -y^r ( E l abcd ®c-a^O ad + E l abcc Vc) 

\c,cZ=3,8,0 c^3,8,0 / 



-crab = - 8 -^§^({M q ,\ a }{M q ,\ b }) 

^0 



^ = ( E T-^e^j^e^ + E 7 a % c M c 2 

;; \c,<2=3,8 c^3,8 

E -M abcd e^e sd + E -^^V (20) 

d=3,8,0 c^3,8,0 / 





Bp 
24F 2 



and, 



1 M 2 
T° M = {[A°,A C ][A*,A'']) 



I abed 



\ E (AtqA^A^A^A^). (21) 



* 2 

perm.{(j} 

Throughout, the terms whose renormalized LECs are set to vanish have not been displayed 
explicitly. It is interesting to note that (p a b does not receive any loop contributions from the 
singlet pseudoscalar mode, implying that F a t, is also free of such contributions. 

It is useful at this point to give the explicit expressions of the masses and decay constants 
at NLO disregarding the chiral logarithms. For the decay constants the above expressions 



lead to 



F n + = F H — (m u + m d ) 

-TO 

-^33 = -F,r+ 

F 88 = F + ° (m M + m d + 4m s ) 

F 00 = F (l + — ) + — (m u + m d + m s ) 
Z d-P 

4L 5j B , s 
-^30 — v^2F 38 



, v 32L 5j B , 
F 80 = — — {m u + m d -2m s 



while for the masses, 



M 2 + = 25 m + 32(2 ^ 2 " L5) ^ 2 

M|+ = 5 (m u + m s ) + 8 ^ 2L8 ~ L Bl(m u + m s ) 5 

^0 

M^o = B {m d + m s ) + %{2L \ l ~ Ls) B 2 {m d + m s ) 2 



M 2 3 = 2B m+ 1 8 b, Bl{ml + m 2 d ) 

^0 

M 8 2 8 = ^ (m + 2m s ) + ^ 2 4|^-^ J Bo 2 K + ^ + 4m 2 ) 

M 2 = M 2 (l - Ai) + mL "~ Lb) Bl{ml + m 2 + m 2 ) 
2 

+ -(1 + p)B (m u + m d + m s ) 
Ml, = -±=B (m u - m d ) + 16(2 ^ " ^ B 2 (ml - m 2 ) 



V3 UV " " V3F 2 



^3 2 o = -\j\{l + ^)B,{m d -m u ) 

16(2L 8 -L 5 )_/2 q2 , , 2 



+ ; q2 — vi^(<-<) 

M 8 2 = ^?(l + |)5 (m u + m d -2m s ) 

1 16(2L 8 -L 5 ) V2 2( 2 2 2 . 

H Y 2 o- B 0\ m u + m d - 2m s)i 



A/f 2 

where p = -A x + 2A 2 - 8L 5 ^ 







The second class of NLO corrections can be grouped into a single term contained in the 
0(p 6 ) odd-intrinsic parity WZ Lagrangian[l2"| — 



= -i7rat 1 ( X _Q 2 )FF 
where x~ = v^X u ^ — ux^u with u = VU (24) 

(There exists a second term |TJ|] that, upon the explicit inclusion of the singlet pseudoscalar 
meson, becomes sub leading in 1/N C and is therefore neglected.) The low energy constant t\ 
has vanishing /5-function and its value can be estimated by means of a QCD sum rule for 
the general three-point function involving the pseudoscalar density and two vector currents 
and saturating the spectral function in the hadronic sector with the states indicated in Fig. 
(2), yielding^, M- 



^-k {F ° + M ] - (25) 

Here the F 2 contribution is determined by the masses and decay constants of the vector 
mesons (the vector meson mass is naturally taken to be my ~ m p ) and is represented by Fig. 
(2b), while the contribution proportional to r is determined by excited pseudoscalars, such 
as the 7r'(1300), and is represented by Fig. (2a). This latter contribution can be estimated 
within a model fl^ , p0| and its expected size is at most one third of the magnitude of the 



vector meson contribution. Since this is similar to the level of uncertainty expected in the 
sum rule evaluation, the r-piece will be disregarded henceforth. As shown by the numerical 
analysis below, the effects on the 7r° width due to the C^J 1 with t\ as estimated above are 
of similar magnitude to the rest of the NLO corrections and in the range of 0.5 %. 
At this point the NLO two-photon amplitudes can be explicitly given: 

(77 I O = -ia {^C a F~i + ^ tl Q- m {{\\ M q }Q 2 )^j (77 I FF | 0) (26) 

and the term proportional to £1 can be obtained in two equivalent ways, either by determining 
the contribution from C^ z to the matrix elements 

(77|g75{A a ,M,}<?|0> 

in Eq. [| or by directly calculating the contribution to (77 | rr a ) due to that effective 

T.acrrfmcnfm TNTntp tl~mt tVip pnntril~intinm from hai tVip cmmp opaline in AT fir thp lppHino - 



(a) 



(b) 



FIG. 2: Hadronic contributions to t\\ tt' denotes excited pseudoscalar mesons and V denotes 
vector mesons. 

one. In Eq. [26] the factor I?o({A a , A4 q }Q 2 ) can be expressed using the LO mass formulas, 
namely 



B Q ({X a ,M q }Q 2 ) 



l(3M 2 + 5(M 2 K+ - M 2 R0 )) 
^(7M 2 + M 2 K+ -hM\ 
lJl(2M* + 2M 2 K+ -M? 



a = 
a = 
a 



(27) 







At NLO the extraction of the ratio R = m s /(m^ — m u ) that characterizes the size of 
isospin breaking should be improved by including NLO corrections to Dashen's theorem. 
Over time several works have shown that the corrections are sizeable. From the early works 
of Donoghue, Holstein and Wyler [21] and Bijnens [22], and more recent works P3[] , it is well 
established that the mass difference M#o — Mk+ left after subtracting the EM contribution 
is larger than the one predicted by Dashen's theorem. While Dashen's theorem predicts 
Mft-o — Mk+ = 5.25 MeV, after the corrections have been implemented it is estimated that 
M K o — Mk+ = 6.97 MeV. The mass difference is slightly smaller if the chiral logarithms 
are neglected following the approach of this work. In such a case M K o — M K + = 6.47 MeV. 
These corrections to Dashen's theorem translate naturally into a smaller value for R, as 
shown in the analysis that follows. There is additional evidence that R is overestimated 
by applying Dashen's theorem, and that comes from the decay r] — > 7r + 7r~7r°. At lowest 
order the decay rate, which is proportional to (m^ — m u ) 2 , is found to be a mere 66 eV 
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At NLO in the chiral expansion the rate is increased to 167 ± 50 eV jnj, while dispersive 



analyses PS] , ^] give 209 ± 56 eV, which is still substantially below the experimental value. 
Clearly one way to make up for the difference is to increase — m u . Because of the large 
uncertainties in both the experimental and theoretical side it is difficult to be precise, but 
it seems that the increment implied by the violations to Dashen's theorem mentioned above 
is in line with the enhancement required to explain the observed rj — > 7r + 7r~7r° width. In 
principle a precise measurement of the tt° width could provide an independent determination 
of R. However, as shown by the results of the present analysis, the 7r° width is affected by 
the corrections to Dashen's theorem only at the level of 0.5%, which is unfortunately well 
below the experimental error of 1.4% aimed at by PRIMEX and about the same as the 0.6% 
uncertainty due to the 0.3% error in the the experimental value of F n + 0. 

There are nine low energy constants to be determined — F , B mi (i = u, d, s), M , L 5 , L 8 , 
Ai, and A2 and these can be found by solving for the observables: F n + = 92.42 ± 0.25 MeV, 
F K+ = 113.0 ± 1.6 MeV, M n o = 134.976 MeV, M v = 547.30 MeV, M n > = 957.78 MeV, 
M K o = 497.78 MeV, M K o - M K + (which as mentioned before, is 5.25 MeV at LO, while 
at NLO and disregarding (including) the chiral logarithms in the corrections to Dashen's 
theorem is 6.47 MeV (6.97 MeV), T^ 77 = 464 ±45 eV, and T v ^ 77 = 4.28 ±0.34 keV. Note 
that the tenth LEC t\ cannot at this stage be extracted phenomenologically and thus its 
value is taken according to the estimate made above. 

In Table I the second column displays the LO results, and the next three columns display 
three different NLO fits, namely: 

i) NLO#l includes terms of 0{p & ) and 0{p A /N c ) in the decay amplitude — i.e. chiral 
logarithms are omitted. 

ii) NLO#2 includes chiral logarithms, which are 0(p 6 /N c ), and the renormalization scale 
fi is set equal to M v . 

iii) NLO#3 is identical to NLO#l but sets h = —i.e. excludes the 0(p 6 ) WZ contri- 
butions. 



TABLE I: Results for the LECs for several different fits, the LO fit and three NLO fits as explained 
in the text. The LECs that depend of the QCD renormalization scale ixqcd correspond to the 
limit ixqcd — ► oo. Although all the LECs are the renormalized ones, only those in the NLO fit #2 
depend on the chiral renormalization fi that is taken to be equal to M v . 



LEC 


LO Fit 


NLO #1 


NLO #2 


NLO #3 


F [MeV] 


92.5±0.6 


90.73 ± 0.47 


84.54 ±0.85 


90.69 ± 0.38 


M [MeV] 


848±40 


1047 ± 5 


1142 ± 33 


1044 ± 4 


2B m [GeV 2 ] 


0.0366 ± 0.0001 


0.03656 ± 0.00001 


0.0362 ± 0.0002 


0.03648 ± 0.00001 


B (m d -m u ) [GeV 2 ] 


0.0235 ± 0.0006 


0.0237 ± 0.0002 


0.0245 ± 0.0006 


0.0244 ± 0.0002 


B Q m s [GeV 2 ] 


0.236 ± 0.006) 


0.2349 ± 0.0006 


0.197 ±0.002 


0.2291 ±0.0005 


2L 5 + L 8 





(5.26 ±0.01) x 10~ 3 


(6.3 ± 1.1) x 10~ 3 


(5.44 ±0.07) x 10~ 3 


2Lg — L5 





(0.8 ±0.9) x 10~ 5 


(-0.53 ±0.04) x 10~ 3 


(0.11 ±0.01) x 10~ 3 


Ai 





0.19 ±0.01 


0.29 ± 0.04 


0.209 ±0.006 


A 2 





0.74 ±0.02 


1.4 ±0.4 


0.81 ±0.02 



It should be noted that the mass of the 7/ in the chiral limit and at NLO in 1/N C is given 
by y/1 — A1M0 ~ 940 MeV, which is slightly high leaving not enough room for the piece 
linear in the quark masses. This linear contribution is suppressed by the rather large value 
of A 2 which leads to a very small value of 1 + p. This cancellation between the leading and 
subleading in 1/N C contributions seems to indicate some difficulty with the 1/N C expansion 
for the masses. Ther first manifestation of this problem is of course in the problem found 
with the 77 and 7/ masses at LO. Although this problem has a minor impact on n° width, it 
certainly deserves further study. It should be noted that the mass difference M v + — M n o has 
been given as input in fit #2 since its value emerges as too large if left unconstrained. The 
reason why it is too large can be traced back to the chiral logarithms generated by the 7/. 
It seems, therefore, that requiring the subleading renormalized LECs to vanish is not such 
a good approximation when such chiral logarithms are included. 

Table II lists the associated predictions for various quantities of interest, in particular the 
7T° width. It is evident from the NLO fits that T^-o^^ is rather stable and always within 
1% of the leading order result. This is within the expected range of the NLO corrections. 
As shown by comparison of the first and third NLO fits, the correction from the 0(p 6 ) WZ 



TABLE II: results implied by the different fits displayed in Table I. The * indicates that the 
quantities are inputs. 





LO Fit 


NLO #1 


NLO #2 


NLO #3 


I>_* 77 [eV] 


8.08 


8.10 


8.16 


8.14 


r?7^77 [eV] 


565 


464 * 


464 * 


464 * 


r^'_» 77 [keV] 


5.1 


4.28 * 


4.28 * 


4.28 * 


M n + - M^o [MeV] 


0.32 


0.24 


0.16 * 


0.21 


m s /rh 


25.9 


25.7 


21.7 


25.1 


R = m s /(m d - m u ) 


45.3 


36.6 


30.9 


37.5 


Os [°] 


1.57 


1.51 


1.88 


1.40 


Os [°] 


-0.56 


-0.68 


-0.94 


-0.59 


Oo [°] 


-18.6 


-10.6 


-8.7 


-12.2 



Lagrangian CtJJ 1 reduces the 7r° width by 0.5%, a magnitude in line with the fact that such 
a correction is controlled by the ratio m u ^/A x . The chiral-logarithm contributions to the 
amplitudes as shown by fit #2, provide an increase of order 0.5% to the tt° width. Since 
these are subleading corrections of 0(p 6 /N c ) to the decay amplitude, they are somewhat 
larger than the 0.2 — 0.3% expected from the ratio m u ^/ (N C A X ) that determines them. This 
problem is similar to the one with the pion mass difference just mentioned; in this case the 
rf loops affect the mixing angles producing a larger than expected correction to the rate. 
Indeed, turning off the chiral logarithms generated by the rf the tt° width is essentially 
identical to the result in fit #1. 

It is important to note that the mixing angles are substantially modified at NLO: in 
particular, the 7T° — r\ mixing angle is found to be e ~ 8 3 + 6 8 = 0.8° — 0.9° in the three fits, 
which is ~ 10 — 20% smaller than the LO result of 1°, but still larger than that obtained at 
LO within SU(3) and given in Eqn. |13[ The latter is chiefly a consequence of the corrections 
to Dashen's theorem. The tt° — 7/ mixing e goes from 0.5° at LO to approximately 0.3° at 
NLO; finally, the 77 — rj' mixing angle is dramatically reduced to about —10° from its LO 
value of —18.6°. In view of these substantial corrections to the mixing angles, the stability 
of the 7T° width is nontrivial: besides the corrections due to the 0(p 6 ) WZ Lagrangian, the 
decay amplitude is determined by the decay constant matrix F aa , which is affected by the 



mixing of states as well as by the NLO corrections contained in the decay constant matrix 
F a b given in Eqn. [T^. It turns out that the entries in F^ 1 — namely F~ Q — affecting the n° 
amplitude remain stable well within the natural size of the NLO corrections. 

In order to asses the theoretical uncertainty of the analysis of the 7r° width, an estimate 
of the magnitude of EM corrections should also be given. Such corrections are of order 
a/27r, which puts them in the 0.2 — 0.3% range. Note also that the value of F n being used 



is that of F n + which has an EM correction. This correction has been studied p7| an d is 
given by 5em-Ft ~ K47raF where the low energy constants that determine the coefficient k 
can be estimated in a resonance saturation model that gives |ac| ~ 10 -2 , thus leading to the 
estimate |<We-^7t| ~ 0.1 MeV which is well within the experimental uncertainty in the value 
of F n + , and implies a correction to the 7r° width consistent with the range mentioned above. 
It seems, therefore, safe to neglect electromagnetic corrections at the level of precision of 
this work. 

Although the rj — rj' complex is not the primary focus of this work, the analysis carried 
out illuminates crucial aspects of this system. At LO the description is rather poor, in 
particular because the two-photon widths depart quite substantially from the experimental 
world averages. Indeed, rf£ = 613 eV versus the world average experimental value of 
464 ± 45 eV, and r^ 77 = 4.86 keV vs 4.28 ± 0.34 keV. These latter disagreements are 
mostly due to the large r\ — rf mixing angle that results from the LO mass formulas, and 
the fact that at LO all decay constants are set to be equal to F n . At NLO the scheme that 
emerges is the one already found in other works |14|], where the mixing angle of the pure 
U(3) states is in the proximity of —10° rather than the —20° obtained in LO, and where the 
decay constant matrix can be parametrized by means of two angles and two decay constants 



fi~4] , |28|| . Following the conventions and notation in[14[], the present analysis gives (quantities 
are denoted in boldface not to be confused with quantities defined heretofore in the text): 
F ~ 116 MeV, F 8 ~ 122 MeV, 9 8 ~ -20° 9 = -2.5° to 0.5°, and the angles 6 and 9 8 
for the three NLO fits are respectively (-0.9°, -19.9°), (2.0°, -19.0°) and (-2.5°, -21.5°) 
. The difference 6q — @s turns out to be between 19° and 21°, to be compared with the 
19° obtained in (in a NLO estimate in that reference a value of 14° is obtained, which 
departs substantially from the one of the current analysis). There exist numerous studies of 
the tj — if complex. It makes sense only to compare results with those using the two-angle 
scheme . Although some of the analyses in these references are purely phenomenological, 



these results are in general in good agreement with the results obtained in this work. 

The quark mass ratios obtained in the different fits deserve comment. The ratio m s /rh is 
in good agreement with the standard value 24.3 ± 1.2 obtained in SU(3) J50]; in fit #2 it is a 
couple of standard deviations smaller, most likely because the chiral logarithms included do 
not represent the full NNLO contributions. In all, this is not surprising as the assumption 
that the low energy constants that are subleading in 1/N C can be disregarded is one of the 
important assumptions in the extraction of the standard ratios. However, a comment is in 
order: it is observed that using LO mass formulas in the NLO results — i.e., expressing quark 
masses in terms of meson masses squared in the NLO terms — leads to a fit that is less 
stable and generates large corrections to the quark mass ratios. The ratio R is smaller here 
than the standard value 42.3 ± 4.5, and this is simply because the corrections to Dashen's 
theorem have been included. The values for R in fits #1 and #3 are about one standard 
deviation smaller than the standard value, while in #2 the chiral logarithms involving the 
rj' loops give a substantial reduction (when these are turned off R is similar to the result in 
the other fits). An interesting observation is that setting m u = leads to an inconsistent fit 
and a value for the 7r° width of 8.5 eV. 

V. CONCLUSIONS 

The decay rate for tt° — > 77 has been calculated within a combined chiral and 1/N C 
expansion. At leading order in the expansion, the isospin-breaking induced mixing of the 
pure U(3) states increases the size of F^o-^ from the value 7.725 eV predicted by the 
lowest order chiral anomaly by more than 4.5%. This effect is largely due to the fact that 
the contributions from mixing with the 7r§ and 7To add constructively, and are of similar 
magnitude. However, at LO the resulting 77, rj' — > 77 widths are found to be too large, and 
in general the fit is quite poor. There is a clear need then for the NLO calculation, both in 
order to improve the results in the 77 — 77' sector and to test the stability of the enhancement 
of F n o_, iy observed at LO. The NLO calculation reveals that the LO result for F^o^^ is 
quite robust, being modified by the NLO corrections by less than 1%. This stability is, 
however, non-trivial. Indeed, as already noted, at NLO the mixing angles are substantially 
affected — the mixing angles e and e are reduced by 10 to 30% (a more dramatic reduction 
of roughly 50% results for the mixing angle 9q). The ir° width, however, is only slightly 



affected because the effects that ultimately determine the corrections to the amplitude are 
in the decay constant matrix F aa shown in Eq. 17. This matrix is affected by the mixing, 
and also receives NLO corrections that reside in F a f,, and apparently the NLO modifications 
to the mixings are partially compensated by the latter corrections in the case of the entries 
relevant to the 7r°. 

The primary source of theoretical uncertainty in the present calculation of r^o^^ resides 
in the value of R, which has an uncertainty of about 15%. Using the empirical formula 
resulting from the results in Table II, T n o^ iy ~ (7.725 + 14.1/ R) eV, the uncertainty in R 
translates into an error in the 7T° width of 0.6%. Other sources of uncertainty are the NNLO 
corrections, of which the chiral-logs are an example, and which should be expected to be 
in the range of 0.2 — 0.3%, and also EM corrections which according to a straightforward 
estimate should be in a similar range. All this together implies an uncertainty of less than 1% 
in the theoretical prediction of r^o^^. Note that there is an overall uncertainty in T^o^^ 
due to the 0.3% error in F n +. The fact that the theoretical prediction for r^o_^. 77 shows the 
4.5% enhancement that can be experimentally observed, and the fact that the experimental 
result with the smallest quoted error [|J lies more than three standard deviations below that 
prediction lend great significance to the upcoming PRIMEX measurement. 

It is evident from the above analysis that no predictions for the r] and rj' two-photon widths 
can be made. Rather, these quantities are inputs, and their precise values do not affect in 
any dramatic way the 7r° width. In a more complete study, wherein the analysis is extended 
to additional processes such as rj — > 7r + 7r~7r°, a more precise experimental knowledge of 
such widths would be necessary. A more extensive analysis would also illuminate the 1/N C 
corrections encoded in the LECs Ai and A 2 and help determine them more precisely. 

As this manuscript was being completed an analysis by Ananthanarayan and Moussallam 



| 3l| was posted where the electromagnetic corrections are studied in detail. Their result for 
the 7T° rate is in agreement with the results in this work. Related work is also being completed 
by Kaiser and Leutwyler [p2 |. 
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